Numerical simulation of boundary-layer transition

The transition to turbulence in boundary layers was investigated by direct numerical solution of the nonlinear, three-dimensional, incompressible Navier-Stokes equations in the half-infinite domain over a flat plate. Periodicity was imposed in the streamwise and spanwise directions. A body force was applied to approximate the effect of a nonparallel mean flow. The numerical method was spectra, based on Fourier series and Jacobi polynomials, and used divergence-free basis functions. Extremely rapid convergence was obtained when solving the linear Orr-Sommerfeld equation. The early nonlinear and three-dimensional stages of transition, in a boundary layer disturbed by a vibrating ribbon, were successfully simulated. Excellent qualitative agreement was observed with either experiments or weakly nonlinear theories. In particular, the breakdown pattern was staggered or nonstaggered depending on the disturbance amplitude

Numerical simulation of separated flows

A new numerical method, based on the Vortex Method, for the simulation of two-dimensional separated flows, was developed and tested on a wide range of gases. The fluid is incompressible and the Reynolds number is high. A rigorous analytical basis for the representation of the Navier-Stokes equation in terms of the vorticity is used. An equation for the control of circulation around each body is included. An inviscid outer flow (computed by the Vortex Method) was coupled with a viscous boundary layer flow (computed by an Eulerian method). This version of the Vortex Method treats bodies of arbitrary shape, and accurately computes the pressure and shear stress at the solid boundary. These two quantities reflect the structure of the boundary layer. Several versions of the method are presented and applied to various problems, most of which have massive separation. Comparison of its results with other results, generally experimental, demonstrates the reliability and the general accuracy of the new method, with little dependence on empirical parameters. Many of the complex features of the flow past a circular cylinder, over a wide range of Reynolds numbers, are correctly reproduced

Numerical simulation of vortex breakdown by the vortex-filament method

The vortex filament method was applied to the simulation of vortex breakdown. The principal vortex region was represented by multiple filaments, and an axial velocity component was induced by a spiral winding of the filaments. First, an accuracy check was performed for a cylindrical swirling flow with simple analytical expressions for the axial and theta velocities. The result suggests that the flow field is simulated to any accuracy by increasing the number of filaments. Second, an axisymmetric type vortex breakdown was simulated, with experimental data serving as upstream conditions. The calculated axial and theta velocity contours show the breakdown of the vortex, including a rapid change in the vortex core, followed axially by a recovery zone and then a second breakdown. When three dimensional initial data are used the second breakdown appears to be of the spiral type in correspondence with experimental observations. The present method is easily used to simulate other types of vortex breakdown or other vortex flows with axial velocity

Self similarity of two point correlations in wall bounded turbulent flows

The structure of turbulence at a height y from a wall is affected by the local mean shear at y, by the direct effect of the wall on the eddies, and by the action of other eddies close to or far from the wall. Some researchers believe that a single one of these mechanisms is dominant, while others believe that these effects have to be considered together. It is important to understand the relative importance of these effects in order to develop closure models, for example for the dissipation or for the Reynolds stress equation, and to understand the eddy structure of cross correlation functions and other measures. The specific objective was to examine the two point correlation, R sub vv, of the normal velocity component v near the wall in a turbulent channel flow and in a turbulent boundary layer. The preliminary results show that even in the inhomogeneous turbulent boundary layer, the two-point correlation function may have self similar forms. The results also show that the effects of shear and of blocking are equally important in the form of correlation functions for spacing normal to the wall. But for spanwise spacing, it was found that the eddy structure is quire different in these near flows. So any theory for turbulent structure must take both these effects into account

A review of quasi-coherent structures in a numerically simulated turbulent boundary layer

Preliminary results of a comprehensive study of the structural aspects of a numerically simulated number turbulent boundary layer are presented. A direct Navier-Stokes simulation of a flat-plate, zero pressure gradient boundary layer at Re0 = 670 was used. Most of the known nonrandom, coherent features of turbulent boundary layers are confirmed in the simulation, and several new aspects of their spatial character are reported. The spatial relationships between many of the various structures are described, forming the basis for a more complete kinematical picture of boundary layer physics than has been previously known. In particular, the importance of vortex structures of various forms to the generation of Reynolds shear stress is investigated

A numerical study of bifurcations in a barotropic shear flow

In the last few years, more and more evidence has emerged suggesting that transition to turbulence may be viewed as a succession of bifurcations to deterministic chaos. Most experimental and numerical observations have been restricted to Rayleigh-Benard convection and Taylor-Couette flow between concentric cylinders. An attempt is made to accurately describe the bifurcation sequence leading to chaos in a 2-D temporal free shear layer on the beta-plane. The beta-plane is a locally Cartesian reduction of the equations describing the dynamicss of a shallow layer of fluid on a rotating spherical planet. It is a valid model for large scale flows of interest in meteorology and oceanography

Understanding transition and turbulence through direct simulations

Direct simulations consist in solving the full Navier-Stokes equations, without any turbulence model, and describing all the detailed features of the flow. Usually the flows are three-dimensional and time-dependent and contain both coarse and fine structures, which makes the numerical task very challenging in terms of both the algorithm and the computational effort. Most of the work until now has involved spectral methods, which are highly accurate but not very flexible in terms of geometry or complex equations. For that reason, future work will also rely on high-order finite-difference or other methods. Direct simulations complement experimental work, and both contribute to the theory and the empirical knowledge of turbulence. Once such a simulation has been shown to be accurate, the flow field is completely known in three dimensions and time, including the pressure, the vorticity and any other quantity. On the other hand, most simulations to date solved the incompressible equations in rather simple geometries, and direct simulations will always be limited to moderate Reynolds numbers. Extensive simulations have been conducted in homogeneous turbulence, channel flows, boundary layers, and mixing layers. Much effort is devoted to addressing flows with compressibility and chemical reactions, and to new geometries such as a backward-facing step

Airloads on bluff bodies, with application to the rotor-induced downloads on tilt-rotor aircraft

The aerodynamic characteristics of airfoils with several flap configurations were studied theoretically and experimentally in environments that simulate a wing immersed in the downwash of a hovering rotor. Special techniques were developed for correcting and validating the wind tunnel data for large blockage effects, and the test results were used to evaluate two modern blockage effects, and the test results were used to evaluate two modern computational aerodynamics codes. The combined computed and measured results show that improved flap and leading-edge configurations can be designed which will achieve large reductions in the downloads of tilt-rotor aircraft, and thereby improve their hover efficiency

Investigation of Phase-Lagged Boundary Conditions for Turbulence Resolving Turbomachinery Simulations

The present work explores whether phase-lagged boundary conditions can be used to perform scale resolving simulations of turbomachines. The phase-lagged approach considered is based on storing the flow signal, both at the pitch-wise boundaries and the rotor-stator interface, as its temporal Fourier coefficients for a finite number of harmonics. The method is implemented in an in-house CFD solver, G3D::Flow, which can perform both URANS and hybrid RANS/LES simulations. In order to evaluate the performance of the phase-lagged boundary condition, a comparison is made with a sliding mesh simulation on a compressor cascade. Furthermore, the possibility of breaking an error feedback loop generated in the sampling process by including multiple blade passages is also investigated. It is found that this approach greatly improves convergence and accuracy of the sampling

Direct Numerical Simulation and Theories of Wall Turbulence with a Range of Pressure Gradients

A new Direct Numerical Simulation (DNS) of Couette-Poiseuille flow at a higher Reynolds number is presented and compared with DNS of other wall-bounded flows. It is analyzed in terms of testing semi-theoretical proposals for universal behavior of the velocity, mixing length, or eddy viscosity in pressure gradients, and in terms of assessing the accuracy of two turbulence models. These models are used in two modes, the traditional one with only a dependence on the wall-normal coordinate y, and a newer one in which a lateral dependence on z is added. For pure Couette flow and the Couette-Poiseuille case considered here, this z-dependence allows some models to generate steady streamwise vortices, which generally improves the agreement with DNS and experiment. On the other hand, it complicates the comparison between DNS and models